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We examine weak measurements of arbitrary observables where the object is prepared in a mixed 
state and on which measurements with imperfect detectors are made. The weak value of an observ- 
able can be expressed as a conditional expectation value over an infinite class of different generalized 
Kirkwood quasi-probability distributions. "Strange" weak values for which the real part exceeds the 
eigenvalue spectrum of the observable can only be found if the Terletsky-Margenau-Hill distribution 
is negative, or, equivalently, if the real part of the weak value of the density operator is negative. 
We find that a classical model of a weak measurement exists whenever the Terletsky-Margenau- 
Hill representation of the observable equals the classical representation of the observable and the 
Terletsky-Margenau-Hill distribution is nonnegative. Strange weak values alone are not sufficient to 
obtain a contradiction with classical models. 

We propose feasible weak measurements of photon number of the radiation field. Negative weak 
values of energy contradicts all classical stochastic models, whereas negative weak values of photon 
number contradict all classical stochastic models where the energy is bounded from below by the 
zero-point energy. We examine coherent states in particular, and find negative weak values with 
probabilities of f6 % for kinetic energy (or squared field quadrature), 8 % for harmonic oscillator 
energy and 50 % for photon number. These experiments are robust against detector inefficiency and 
thermal noise. 

PACS numbers: 03.65.Ta, 42.50.Dv 



I. INTRODUCTION 

Physics is an endeavor to construct a mathematical 
model of nature. The basic mathematical model of clas- 
sical physics is one in which all dynamical variables have 
definite values simultaneously. With the advent of sta- 
tistical mechanics, a probabilistic description was intro- 
duced in which each dynamical variable would have some 
value with a certain probability. This is the most general 
description that classical physics can provide. A great 
variety of phenomena can be explained in terms of such 
a model. For example, most optical phenomena can be 
described in terms of a model of complex, stochastic am- 
plitudes. It was not until 1977 that this model was found 
to break down in an experiment which demonstrated the 
phenomenon of antibunching [1]. 

In this paper we examine a relatively new type of mea- 
surement known as weak measurements [2]. Weak mea- 
surements may be performed in exactly the same way as 
standard von Neumann measurements, but with a weak- 
ened interaction [3]. In a weak measurement, the pointer 
will on average register the expectation value of the ob- 
servable that is measured. However, when the weak mea- 
surement is conditioned on a second postselection mea- 
surement, the pointer registers the real part of what is 
known as the "weak value" of the observable. These weak 
values have caused a lot of controversy, in particular be- 
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cause they may exceed the eigenvalue spectrum of the 
observable. 

The main objective of this paper is to discuss the con- 
ditions under which a classical description of weak mea- 
surements is possible, taking into account practical ex- 
perimental limitations and possibilities. In particular, 
our purpose is to investigate in further detail the failure 
of providing a classical description of coherent states in 
such experiments, and to propose a feasible experiment 
demonstrating nonclassical properties of coherent states. 

In Sec. II we give a brief review of the Glauber classi- 
cality criterion, and discuss it's possible limitations. 

In Sec. Ill and appendices A and B we consider a gen- 
eral class of weak measurements where the object may 
be prepared in a mixed state, and where the pointer may 
be prepared in an arbitrary mixed state of vanishing cur- 
rent density. We consider detectors of arbitrary quantum 
efficiency, and derive a generalized weak value. 

In Sec. IV we demonstrate that the weak value can 
be expressed as a conditional expectation over an infi- 
nite set of different quasi-probability distributions. These 
distributions can be regarded as generalizations of the 
complex Kirkwood distribution [4] or the standard or- 
dered distribution [5]. We denote such distributions by 
S'-distributions. The Terletsky-Margenau-Hill distribu- 
tion [6, 7], or T-distribution, is the real parts of the S- 
distribution. We find that strange weak values (for which 
the real part exceeds the eigenvalue spectrum of the ob- 
servable) can only exist if the T-distribution takes nega- 
tive values. We demonstrate that a classical model of a 
weak measurement exists whenever the representation of 
the observable coincides with the classical representation 
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of the observable and the T-distribution is nonnegative. 
Wc demonstrate in particular that negative weak values 
of energy contradict all classical stochastic models, and 
that negative weak values of photon number contradict a 
stochastic model where the energy is bounded from below 
by the zero-point energy. 

In Sec. V we propose two feasible practical realizations 
of weak measurements of photon number and energy. 

In Sec. VI we consider in particular weak measure- 
ments on coherent states. We demonstrate that coherent 
states display negative weak values with probabilities of 
16 % for kinetic energy (or squared field quadrature), 8 
% for harmonic oscillator energy and 50 % for photon 
number. We find that these effects are robust against 
detector inefficiency and thermal noise. 



II. THE GLAUBER CLASSICALITY 
CRITERION 

Glauber and Sudarshan have demonstrated that any 
density operator can be expanded diagonally in terms of 
coherent states [8, 9] 



d 2 aP(a) | a) (a 
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The weight function P{a) is known as the P-distribution. 
Furthermore, one may write the expectation value of any 

normal ordered operator Af 0{a),a) as an integral of 

the form 



(TV 6{a\a) )= j d 2 aO{a*,a)P(a), 



(2) 



where a* is the c-number representation of and a is 
the c-number representation of a. Therefore, if the P- 
distribution has the properties of a valid probability dis- 
tribution, one may say that a classical representation ex- 
ists for any normal ordered operator product. This is 
known as the "optical equivalence theorem" [9] . 

The optical equivalence theorem is the basis of the 
Glauber classicality criterion, according to which all 
states for which the P-distribution is a probability distri- 
bution are regarded as essentially classical. If P fails to 
be a probability distribution, the state is considered as 
nonclassical [8, 10-13]. The Glauber criterion is widely 
accepted as giving the most general distinction between 
quantum and classical states. It is the basis of various 
measures of "nonclassicality" . For example, Hillery has 
defined a measure of nonclassicality based on the distance 
in Hilbcrt space between the object state and coherent 
states [14]. Lee has defined a nonclassical depth defined 
as the minimum average number of thermal photons that 
must be added to render the P-distribution nonnegative 
[15]. More recently, Vogel et.al have derived a hierar- 
chy of observable conditions to test the Glauber criterion 
[16, 17]. 



The Glauber criterion must be considered a conjec- 
ture rather than a proven theorem. It is based on some 
plausible arguments. Firstly, due to the multitude of ar- 
guments in favor of coherent states as the only classical- 
like pure states, it is natural to assume that also classi- 
cal mixtures of coherent states are classical. This is in 
fact equivalent to Glauber's classicality criterion, as can 
be seen from Eq. (1). If the P-distribution is a proba- 
bility distribution, the density matrix can be expressed 
as a classical mixture of coherent states. Secondly, a 
nonnegative P-distribution ensures that the whole range 
of s-ordered Wigner distributions are also nonnegative. 
Therefore, the P-distribution is most "sensitive" to non- 
classicality among all s-ordered distributions. 

It is puzzling that a definition of nonclassicality de- 
pends on the ability to represent expectations of normally 
ordered operator expressions as classical expectation val- 
ues over a probability distribution. This is sometimes 
attributed to the fact that normal ordering of operators 
is closely associated to the theory of photo-detection. 

However, we may devise experimental procedures re- 
lated to operator orderings different from normal order- 
ing of annihilation and creation operators, that may dis- 
play clear contradictions between the classical and quan- 
tum descriptions for the same experiment. This is actu- 
ally the case of weak measurements, whose statistics are 
related to T-distributions instead of to the more stan- 
dard s-ordered distributions. We show below that weak 
measurements lead to strange outputs provided that the 
T-distribution takes negative values. Moreover, we will 
show in this paper that a classical stochastic model may 
fail to describe weak measurements even when the P- 
distribution is a nonnegative probability distribution. 



III. WEAK MEASUREMENTS 

A "measurement" comes about when an auxiliary sys- 
tem interacts with an object. By examining the prop- 
erties of the auxiliary system after the interaction, it 
may be possible to extract information about the ob- 
ject. The auxiliary system is frequently called a "mea- 
surement apparatus" or a "pointer system". The basic 
theory of quantum measurement was examined by von 
Neumann in his seminal work on the mathematical foun- 
dations of quantum mechanics [18]. In this work, von 
Neumann represented the interaction between the object 
and the pointer by an interaction Hamiltonian of the form 
(throughout we will use units in which h = 1) 



eS(t) z> ® P. 



(3) 



A short explanation of the terms involved is in place. The 
constant e represents the interaction strength. The inter- 
action is of short duration, represented in idealized form 
by the 5-term. The hermitian observable v belongs to 
the object Hilbert system H s , and is the observable that 
we want to "measure" . The observable P is the pointer 
momentum, and belongs to the pointer Hilbcrt space 7i a . 
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Although seemingly artificially constructed, this interac- 
tion model has served as an archetype of the interaction 
mechanism in quantum measurements. It has been found 
that the conclusions that can be drawn from this model 
are generic to a number of other interaction mechanisms 
(for a closer discussion of the specific properties of this 
interaction Hamiltonian, see Ref. [19]). 

In a standard, projective measurement, the pointer po- 
sition Q, with [Q, P] = i, displays one of the eigenvalues 
of the object observable v after the measurement inter- 
action. It was demonstrated by von Neumann that in 
order to accomplish this, the state of the pointer prior to 
the interaction should have a small position spread [18]. 
The same effect can be accomplished by allowing the in- 
teraction strength e to be sufficiently large [19]. For this 
reason, this type of measurement is frequently called a 
"strong measurement" . 

Until recently, it was thought that strong measure- 
ments are the only useful type of measurements in quan- 
tum mechanics. However, in 1988 Aharonov et.al. pro- 
posed a new type of measurement that they called "weak 
measurements" [2] . Such measurements also employ the 
von Neumann interaction mechanism (3), the difference 
being that the pointer is assumed to be in an initial state 
of large position uncertainty. More specifically, they as- 
sumed that the initial state of the pointer was a Gaussian 
with large spread. 

Recently, it was shown that weak measurements can 
be performed also when the pointer is in an arbitrary 
mixed state, provided that the interaction strength e is 
sufficiently small and the current density of the pointer 
vanishes [3]. This description can be generalized further 
by taking into account finite efficiency of the detectors. A 
detector of finite efficiency can be represented by Positive 
Operator Valued Measure (POVM) (see Appendix A). In 
this case, the weak value of the observable v conditioned 
on an imperfect postselection of the observable <j> on the 
object is (see Appendix B) 



M0) 



Tr(ll0Z>p s ) 
Tr(f[ 0/ S s ) ; 



(4) 



where II^ is a diagonal POVM representing the imperfect 
postselection and p s is the initial state of the object (see 
Eq. (B5)). 



IV. ON CLASSICAL MODELS OF WEAK 
MEASUREMENTS 

Under what circumstances is it possible to find a clas- 
sical representation of a weak measurement? Or put dif- 
ferently, under what circumstances can the outcome of 
a weak measurement be modelled in terms of a clas- 
sical, stochastic model? By answering this question, 
we will also understand under what circumstances weak 
measurements demonstrate nonclassical properties of the 
quantum state under consideration. 



We begin by inserting the definition (B5) for the diag- 
onal postselection POVM into Eq. (4). We then have 

, _ T±UjW IW) |^)(^ \j>p a ] 
" wW ~ Tr[f#'IW)l0')(0'l/3 s ] ' {) 

By performing the trace over any complete set of states 
we obtain 

VwW ~ f#'ii^')(0'l/M0'> ■ {b) 

We employ an arbitrary complete set of states | £), 



/ 



deio«i=i- 



(7) 



We may then write Eq. (6) in the form 
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where 



(010 ' 



(9) 



is a c-number representation of the observable v. In fact, 
it is the weak value for the observable preselected in the 
state | £) and postselected in the state | 0). Also, 



5(0,0 = (010(^1 



(10) 



is a generalization of the Kirkwood distribution for arbi- 
trary observables £ and <j> [4] . The Kirkwood distribution 
is also known as the anti-standard ordered distribution 
[5]. S is also a generalization of the standard ordered 
distribution, which is the complex conjugate of the Kirk- 
wood distribution [5]. In this paper, we will simply de- 
note it by the S-distribution, and we will refer to So as 
the S'-representation of the observable v. 

The ^-distributions are in general complex, and as 
such are quasi-probability distributions. Nevertheless, 
they possess some of the properties of classical joint dis- 
tributions. For example, assuming that both eigenstates 
| (f>) and | £) constitute complete sets, it is straightfor- 
ward to show that they yield correct marginal distribu- 
tions when integrated over either variable, 

(01/5,1^) - J de (ii) 

(£|p s |0 = / d4>S(4>,0- (12) 



It is straightforward to show that also the complex conju- 
gate distribution, S*(<p,£), fulfills such marginality con- 
ditions. 

The .^-distribution can be expressed in the form 



5(0,0=1 (010 
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(00) 



(13) 
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In this form, the ^-distribution is a product of a non- 
negative probability distribution and the weak value of 
the density operator. If | £} or | (f>) are eigenstates of the 
density operator, S will be real and nonnegative. 

Classically, the weak value of an observable is the con- 
ditional expectation of that observable [3]. The expres- 
sion (8) demonstrates that there exists an infinite set of 
representations under which the weak value can be ex- 
pressed as a conditional expectation of a c-number vari- 
able So over an 5-distribution. For each choice of ba- 
sis | £), a different 5-distribution is obtained. However, 
some of these representations bear little resemblance to 
any classical model. For example, in some models the 
S'-representation of the hermitian observable v is com- 
plex. As a basic requirement on a classical model, we 
shall in the following restrict the attention to the subset 
of representations for which hermitian observables have 
real representations, 



Im$»(0,O=O. 



(14) 



We also introduce the T-representation of v and the T- 
distribution [6, 7], 



T(<f>,0 



ReS(0,O- 



(15) 
(16) 



It is straightforward to show that also these provide cor- 
rect marginal distributions, 

(4>\ps\4>) = J d£T(<t>,£), (17) 

(£ \h |0 - J #T(0,O- (18) 

Because of the classicality condition (14), the S- and T- 
representations of the observable are the same, 



mo 

(010 ■ 



(19) 



The distribution of the postselection observable <f>, taking 
into account the finite detector efficiency represented by 
n^, can be found both from the complex S-distribution 
and the real T-distribution through the integrals 

Tt(U^ s ) = JdSdtflliWW,® 



In a weak measurement, the real part of the weak value 
v w is registered by the pointer. Under the assumption 
(14), we may write 

This expression, which reflects the expectation of the 
pointer displacement in a weak measurement, bears for- 
mal resemblance to a classical conditional expectation. 



The difference is that the T-distribution may take nega- 
tive values, and that the T-representation of the observ- 
able may differ strongly from the representation of the 
observable in classical theory. A classical representation 
of the pointer displacement can be said to exist if T,j 
equals the classical representation of the observable and 
if also the T-distribution is nonnegative. 

What determines the sign of the T-distribution? From 
Eqs. (13) and (16) it follows that 



T(<f>,£)=\ (010 



2 D _/ (Z\p.\4>) 



(22) 



Thus, the sign of the T-distribution equals the sign of 
the real part of the weak value of the density operator. 
Therefore, a requirement for nonclassicality is that the 
real part of the weak value of the density operator should 
be negative. 

In the following, we consider two representations that 
may provide classical-like models of weak measurements. 
The first representation, which we shall call the eigen- 
value representation, is found when the basis | is cho- 
sen as eigenstates of the observable v. In this representa- 
tion the S'-and T-representations of the observable v are 
the eigenvalues v, 

Sp((f>,v)=T p (<l>,v) = v, (23) 

By using this representation, we have 



v w {4>) = J duvSj,{u | 0), 
Rev w (<j)) = J dv v T v (v \ 4>) , 



where 



S v (v | 4>) 
T v (v | <f>) 



S v (v, <t>) 

j dvs v (v, <t>y 

T v (v, 4>) 



arc "effective", conditional distributions, and where 



S n {<j),v) 



y iw)s(0» 



(24) 
(25) 

(26) 
(27) 

(28) 



T v {<P,u) = J dcj>> n^'W' » (29) 

are "effective" S- and T-distributions for v and <j>. 

It follows straightforwardly from Eq. (25) that 
"strange weak values" where the real part of v w exceeds 
the eigenvalue spectrum is only possible for quantum 
states for which the T-distribution in the eigenvalue rep- 
resentation takes negative values. But it is of course well 
known that classical models may allow observables to 
exceed the eigenvalue spectrum of the observable. This 
means also that "strange" weak values may sometimes by 
supported by a classical model. In particular, this may 
be the case for observables with a discrete spectrum. 
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The other c-number representation that we shall con- 
sider is the phase space representation. In this case, the 
basis | £ ) should be chosen as eigenstates of the observ- 
able canonically conjugate to the postselection observable 
<j>. This may provide the possibility of a comparison with 
a classical phase space description of weak measurements. 

We give a couple of examples illustrating the use of 
these two representations. First, we consider a weak mea- 
surement of a field squared quadrature (or kinetic energy 
for material particles) v — p 2 postselected on the canoni- 
cally conjugate field quadrature (or position for material 
particles) (f> = q [20, 21]. In this case, the eigenvalue rep- 
resentation and the phase space representation arc one 
and the same, since q and p are canonically conjugate 
variables. From Eq. (19) follows that both the eigen- 
value representation and the phase space representation 
of p 2 postselected on q is 

T p (q,p)=p 2 . (30) 

Obviously, Tp2(g,p) > 0. From Eq. (25) and the classi- 
cality assumption T(q,p) > follows the inequality 

Re(p 2 J(q)>0. (31) 

In this case, a "strange" negative weak value implies fail- 
ure of the classical model where p 2 takes the positive 
continuum. It can be noted that a negative weak value 
of kinetic energy p 2 contradicts all stochastic c-number 
models where kinetic energy takes only nonnegative val- 
ues, even models where the T-representation of kinetic 
energy may differ from the expression (30). 

Consider next a weak measurement of the energy of a 
harmonic oscillator (we use units so that u = 1) 

H=\{p 2 + q 2 )=n+\, (32) 

assuming postselection on one of the quadratures (or po- 
sition for a material particle). In the eigenvalue repre- 
sentation 

T fl (q,n) = n+\. (33) 

Obviously, T^(q,n) > \- In the eigenvalue representa- 
tion, the T-representation of the Hamiltonian is bounded 
from below by the zero-point energy. From Eq. (25) and 
by imposing the classicality criterion T(q, n) > 0, we may 
derive the inequality 

Re(H w )(q) > 1 (34) 

Violation of this inequality implies "strange" weak val- 
ues. It can only take place if the T-distribution T(q, n) 
takes negative values. We can split the strange val- 
ues into two categories, (i) Outputs of the form i > 
Rc(H w )(q) > contradict exclusively those classical 
stochastic models for which the excitation number takes 



positive values n > 0, or, equivalently, for which the en- 
ergy is larger than or equal to 1/2. These outputs do not 
contradict classical models where the energy takes pos- 
itive values starting from zero, (ii) On the other hand, 
outputs of the form > Re(H w )(q) contradict all the 
classical stochastic models where the energy takes posi- 
tive values. 

In the phase-space representation, a weak measure- 
ment of the energy of a harmonic oscillator postselected 
on one of the quadratures calls for the T-representation 

T k (q,p) = \(v 2 + a 2 )- (35) 

This representation of harmonic oscillator energy coin- 
cides with the classical representation, and therefore per- 
mits an investigation of the limits of classical models. 
Obviously we have T^(q,p) > 0. By assuming that also 
the state representation is classical, T(q,p) > 0, one may 
derive from Eq. (21) the inequality 

Re(H w )(q) > 0. (36) 

This inequality should be compared with inequality (34). 
We may say that the phase-space representation high- 
lights just the second category (ii) above: Violation of 
inequality (36) rules out all the classical stochastic mod- 
els where the energy takes positive values. Note that the 
probability of infringing inequality (36) and the proba- 
bility of the second category (ii) above are equal, since 
the probability of Rc(H w )(q) < does not depend on 
the representation. On the other hand, the probabil- 
ity that Re(H w )(q) < is less than the probability of 
Re(H w )(q) < 1/2 in agreement with the fact that the 
first one excludes a larger class of classical models. 

We have not mentioned the correspondence principle 
here. Our purpose is to investigate under what condi- 
tions quantum mechanics can be reproduced by a classi- 
cal stochastic theory. This is not related directly to the 
classical limit of quantum mechanics. 

V. WEAK MEASUREMENT OF PHOTON 
NUMBER 

The realization of weak measurements requires the 
coupling of the system to be observed with auxiliary de- 
grees of freedom. The output of the weak measurement 
is inferred from measurements carried out on the auxil- 
iary system and on the object system itself. In this work 
we will consider the weak measurement of two observ- 
ables with nonnegative spectra. The boundedness of the 
spectra is mandatory in order to reveal the appearance 
of strange values. These observables are v — p 2 and the 
number operator v = a^a, where a = (q + ip)/V2 and 
q, p, with [q,p] — i are the standard position and linear 
momentum or field quadratures. For both examples of v 
we will consider the same postselection strategy given by 
the measurement of the operator (f> = q. 
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b. nonlinear crystal 



v 
■0 



quadrature 



quadrature 

FIG. 1: Scheme for a weak measurement via cross-Kerr cou- 
pling in a nonlinear crystal. 



In this section we propose two simple and feasible 
schemes for the weak measurement of a*d conditioned 
on the measurement of q in the field of quantum optics, 
where a^a represents the number of photons and q is a 
field quadrature. The two possibilities involve different 
realizations of the auxiliary system. These are another 
field mode and two-level atoms. 



A. Coupling to a field mode 

Let us assume that the auxiliary system is another field 
mode with complex amplitude operator b. A suitable 
coupling between the system a and the auxiliary variables 
allowing a weak measurement of the number operator a' a 
is of the form 



(37) 



This coupling can be achieved in practice by propaga- 
tion of both field modes in crystals with nonlinear optical 
properties (cross-Kerr interaction) so that e is propor- 
tional to the nonlinear susceptibility of the medium and 
the length of the crystal. This coupling causes a phase 
shift of the mode b proportional to the photon number 
in mode a, that can be then detected simply by mea- 
suring a quadrature of the mode b (homodyne detection) 
[22, 23]. This example has the advantage that nonlinear 
effects are usually very weak, so that the requirement 
e — * for a weak measurement is naturally satisfied. In 
Fig. 1 we outline the scheme of the weak measurement of 
photon number via cross-Kerr coupling conditioned on a 
quadrature measurement. 



B. Coupling to a two-level atom 

The weak measurement of the photon number can also 
be carried out by coupling the field mode to a two-level 
atom with internal energy levels |±). If the frequency 
of the field and the resonant frequency of the atom are 
detuned enough the atom-field interaction Hamiltonian 
becomes [22, 23] 



H, 



ea)a<j z 



(38) 



where a z = |+) (+| — |— ) (— |. In this case we have that e is 
inversely proportional to the detuning, so the condition 
e — > can be easily achieved. 



The atom-field interaction causes a phase shift of the 
coefficients of an atomic superposition of the states |±). 
The shift is proportional to the photon number and can 
be detected by measuring, for example, the observable 



^ = !+)<-! + !+)(-!• 



(39) 



The measurement of a x can be carried out by determin- 
ing the population of the levels |±) after applying to the 
atom a resonant pulse transforming the eigenstates \±) x 
of a x into |±). 

The condition of null current density (B9) becomes in 
this case 



x (±\ (<7 z p a + PaCT z ) |±)x = 0, 

which is verified always provided that 



z 'y u y) > 



where s x , s y are real constants such that s x 



(40) 



(41) 



Sy < I and 



<rv = i(\+)(-\-\+)H)- 



(42) 



Besides quantum optics, this scheme can be also im- 
plemented in the context of trapped ions, where a would 
represent the complex amplitude of the one-dimensional 
harmonic motion of the center of mass of the trapped 
ion, and |±) are two internal levels of the same ion. 



VI. COHERENT STATES 

In this section, we turn to the study of weak mea- 
surements on coherent states, taking into consideration 
the effect of thermal noise and finite detector efficiency. 
It is of particular interest to study possible nonclassical 
properties of coherent states, since coherent states are 
the only pure states that satisfy the Glauber criterion of 
possessing a nonnegative P-distribution. 

We shall consider the thermalized coherent state (also 
known as the displaced thermal state) 



= D(a)fHh&(a). 



(43) 



Here plh is the density operator for the thermal state 
and D{a) is the displacement operator. This state has a 
nonnegative P-distribution [24] 



1 



imth 



-17- 



/n t h 



(44) 



where a is the coherent amplitude when nth vanishes and 
nth is the expected thermal photon number when a van- 
ishes. Since the P-distribution is nonnegative, this state 
is essentially classical according to the Glauber criterion. 

In this section, we consider weak measurements with 
postselection on position. The phase space for this exper- 
iment consists of position and momentum. In the phase 
space representation, the ^-distribution for this state is 



S(q,p) = S t h(q -a r ,p- oj), 



(45) 
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where a = (a r 



10L 



Sth{q,p) 



exp 



2< h (p 2 + g 2 )-2ip g 



TT^l + 4(7. 



is the ^-distribution of a thermal state [25] and 



°th = n th + g 



(46) 



(47) 



is the variance of each quadrature for the thermal dis- 
tribution. Clearly, the T-distribution for this state takes 
negative values. It is worth emphasizing that the lack of 
positivity persists for every n t h, in sharp contrast to the 
case of s-ordered distributions for which there is always 
a value of n t h that renders the distribution positive. 

In the following, we assume an imperfect measurement 
of the operator q represented by a Gaussian postselection 
POVM (see also Eq. (B5)) 



IW) 



1 



27T(7„ 



(48) 



where the width a n is determined by detector efficiencies. 
For example, for a homodyne detector, a 2 = (1 — rj)/ (2rj), 
where 77 is the quantum efficiency of a single detector. We 
may then derive an effective marginal distribution for the 
postselection observable 



p v (q) = Tr 



(49) 



A. Negative weak value of kinetic energy 

In this subsection, we consider weak measurements of 
the observable p 2 , which is essentially the kinetic energy 
[20] . The weak value of p 2 conditioned on the measure- 
ment of q was found to take negative values for coherent 
states in Rcf. [21]. Here, this treatment will be general- 
ized to include the effects of thermal noise and imperfect 
detectors. 

Since we are considering weak measurements of p 2 , the 
phase space representation and the eigenvalue represen- 
tation of the 5-distribution coincide. It is then useful to 
define an effective ^-distribution by using Eq. (28) and 
(46). The result is 



SJq,p) = Sf{q 



,P-OLi), 



with 



Si h (q,p) 



exp 



-2ipq 



(50) 



(51) 




FIG. 2: The probability for observing a negative weak value 
of p 2 as a function of the detector efficiency r\ and the average 
thermal occupation number nth- It is assumed that at = 
0, whereas a r may take any value. This effect therefore is 
present even for macroscopic excitations a T . 



being the effective S'-distribution. Using Eq. (25) and 
(51) we find that 

Re i = l + 4(a t 2 + a t 2 JK\ + ^) _ (g - a r ) 2 



The roots of this polynomial are 



(52) 



q± = ^r±^/(a 2 h + a 2 )[l + 4(af + „ 2 h )(a 2 h + a 2 )}. (53) 

The probability of postselecting a position q which on 
average gives a negative weak value of p 2 then is 



P [Re(p 2 ) w < 0] 
The result is 



1 



9+ 



dq p v {q). 



(54) 



P [Re(p 2 ) w < 0] =erfc^/- 



2(« 2 



' th 



(55) 

The complementary error function erfc(ir) is monotoni- 
cally decreasing. Therefore, this probability is maximized 
when a 2 , a 2 h and a 2 are as small as possible. The state 
can be chosen so that a L — 0. For a perfect detector 
(Tjj = 0. The quadrature variance a 2 h is bounded below 
by 1/2. Therefore, the maximum probability for a nega- 
tive weak value of p 2 is erfc 1 w 0.16 [21]. It is particularly 
interesting to note that the effect is independent of the 
real part of the amplitude, a r . A negative weak value 
therefore can be observed for a macroscopic occupation 
of the mode. 

This probability has been plotted as a function of de- 
tector efficiency r\ and thermal occupation number nth in 
Fig. 2. As noted in section IV, the negativity of Re(f>^) 
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contradicts all classical stochastic models. A weak value 
of kinetic energy might be observed for material parti- 
cles at low temperatures. However, a realizable quantum 
optical experiment is not known to the authors. To in- 
vestigate a feasible quantum optical experiment, we turn 
to the weak measurement of photon number and energy. 




B. Negative weak value of energy 

In this subsection we study weak measurements of en- 
ergy conditioned on the postselection of a quadrature ob- 
servable. As demonstrated in Sec. IV, a negative weak 
value of energy for a harmonic oscillator contradicts clas- 
sical stochastic models. 

By combining Eqs. (21), (46) and (48) we find that 



Re[H w (q)) = aq 2 + bq + c 



where 



4< 



1 



8(^+^)2' 

"r(4<T t Vg + 1) 

4(4* 



'/,/), 



a[ 1_ 

2 + 8( 



^(4< - 1) 



(56) 

(57) 
(58) 
(59) 



'th 



Two real roots q± exist provided that b 2 > 4ac (only 
one if a = 0). This establishes a necessary condition to 
be fulfilled by a t h, cr n , a r , and on for the existence of 
negative values for Re [H w (q)] . If this is satisfied, the 
probability of observing a negative Re (H w ) is 



P[Re(H w ) < 0] 



dqp v (q). 



(60) 



For ideal detectors a n = and vanishing thermal noise 
nth — this probability can be written as 



1 



P[Re(H w ) < 0] = - erfc 



1 



- 07 



(61) 



This probability has been plotted in Fig. 3. It reaches a 
maximum at a r = 1 and on = 0, at which the probability 
is (1/2) erfc 1 w 0.08. This is half of the probability for 
observing a negative weak value of kinetic energy. 

A more complex expression for the probability can be 
obtained for arbitrary detector efficiency and thermal ex- 
citation. This has been plotted in Fig. 4. We see that a r 
now must reach a minimum value in order to see nonclas- 
sical negative weak values. In Fig. 5, the probability is 
plotted as a function of detector efficiency rj and average 
thermal occupation number nth- 

Negativity of weak values of kinetic energy and of total 
energy both contradict a classical, stochastic model of 
light. Negative weak values of kinetic energy persists 
also for macroscopic coherent amplitudes. This does not 
occur for the weak value of energy. On the other hand, 
a weak measurement of energy has a feasible practical 
measuring schemes as outlined in Sec. V. 




FIG. 3: The probability for observing a negative weak value of 
energy as a function of the coherent state amplitudes ot r and 
ext. We assume perfect detectors a v = and no thermal noise 
nth — 0. The probability has a maximum at a r = 1, ctj = 0, 
at which the probability is (1/2) erfc 1 w 0.079. 




FIG. 4: The probability for observing a negative weak value of 
energy as a function of the coherent state amplitudes a r and 
«i. We assume detectors with quantum efficiency 77 = 0.7 and 
thermal occupation number n t h = 0.3. There is a minimum 
a r required in order to see nonclassical behavior. 



C. Negative weak value of photon number 

It was demonstrated in Sec. IV that a negative weak 
value of photon number contradicts a classical stochastic 
model where the energy is bounded from below by the 
zero-point energy. We study this further here. 

The weak value of photon number is simply 



n w (q) = H w (q) - -. 



(62) 



Two real roots q± of n w exist provided that b 2 > 4a(c - 



2.5 



FIG. 5: The probability for observing a negative weak value of 
energy as a function of the detector efficiency r\ and the aver- 
age thermal occupation number nth- The coherent excitation 
is assumed to be a r = 1, a, = 0. 

1/2). In the case that this condition is satisfied the prob- 
ability of observing a negative Re is 

P[Re(n w ) < 0] = T + dq Pri (q). (63) 

Jq_ 

For ideal detectors and vanishing thermal noise this prob- 
ability can be written as 

1 fa 2 + a 2 \ 
P [Re(n w ) < 0] = - erfc {jT^j) ' (64) 

This function has been plotted in Fig. 6. The probability 
is always maximized by letting on — > 0. It approaches a 
maximum of 0.5 for vanishing a r . However, it has a 
singularity in a r = 0, and actually vanishes in this point. 
Thus, there is zero probability of observing a negative 
weak value of h for the vacuum state. There must be a 
finite small coherent amplitude to see this. 

A more complex expression is obtained for finite detec- 
tor efficiency and finite thermal noise. The probability 
has been plotted in this case in Fig. 7. We see that a r 
now must reach a minimum value in order to see non- 
classical negative weak values. 

The probability has been plotted as a function of de- 
tector efficiency n and the average thermal occupation 
number nth in Fig- 8. We see that for finite thermal 
occupation number nth there is a lower bound on the 
detector efficiency n to see nonclassical behavior. This 
bound vanishes when n t h — > 0. 

VII. CONCLUSIONS 

In this work we have examined a very general form of 
weak measurements focusing on the emergence of non- 
classical features. We have shown that the appearance 



FIG. 6: The probability for observing a negative weak value 
of the photon number n as a function of the coherent state 
amplitudes a r and at. We assume perfect detectors a v = 
and no thermal noise nth = 0. The probability increases with 
decreasing [ a |, but has a singularity at | a |= 0. 




FIG. 7: The probability for observing a negative weak value 
of the photon number n as a function of the coherent state 
amplitudes a r and Qi. We assume detectors with quantum 
efficiency 77 = 0.7 and thermal occupation number nth = 0.3. 
We see that there is a minimum a r required in order to see 
nonclassical behavior. 



of strange weak values is equivalent to the existence of 
negative values for a generalized Terletsky-Margenau-Hill 
distribution. 

We have presented some feasible practical implemen- 
tations of this kind of measurement in the field of quan- 
tum optics focusing on the weak measurement of photon 
number and energy. We have demonstrated that nega- 
tive weak values of energy contradict all classical models 
of light, and that negative weak values of photon number 
contradict a classical model where the energy is bounded 
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FIG. 8: The probability for observing a negative weak value of 
the photon number n as a function of the detector efficiency rj 
and the average thermal occupation number nth- The coher- 
ent excitation is assumed to be a r — 0.1, at — 0. A minimum 
detector efficiency (and a maximum thermal occupation num- 
ber) is required to observe a negative weak value. 



from below by the zero-point energy. 

As a particular but striking enough example we have 
considered weak measurements on coherent states. We 
have found that negative weak values can be observed 
with a probability of 16 % for kinetic energy (or squared 
field quadrature), 8 % for harmonic oscillator energy and 
50 % for photon number. 

We have analyzed the persistence of the effect under 
practical experimental conditions by considering degrad- 
ing imperfections such as the presence of thermal fluc- 
tuations and the use of inefficient detectors. All these 
results confirm the possibility of a practical observation 
of nonclassical effects for states previously considered as 
firm examples of classical behavior. 



APPENDIX A: IMPERFECT DETECTORS 



The POVM should provide a resolution of the identity 
operator 



n* = i. 



This implies that 



n rf 



i. 



(A3) 



(A4) 



Thus, IL,f,(<p') should be a normalized distribution over </>. 
The probability distribution for the observable 0, taking 
into account the imperfect detector represented by II^, 
is 



Ti- 



ll,/ 



dfill^'M \p\<f>') 



(A5) 



We require that the imperfect detector should give the 
same expected reading as a perfect detector. This condi- 
tion of unbiascdncss can be written as 



y"d0 0Tr [ft^p] = J d<f><t>((j) | 
and it implies that 

4>'= I ddpdp n (0'). 



p\<P), (A6) 



(A7) 



Thus, unbiasedness is equivalent to requiring that the 
parameter <fi' should be the expectation value of the dis- 
tribution . 

A typical diagonal POVM satisfying the properties 
(A4) and (A7) is the Gaussian 



i 



= -(0-0') 2 /(2t 2 ) 



V27 



(A8) 



For this POVM, and for the subset of POVM's where 
H^,((j)') is a function of (0 — (f>') only, the transformation 
(A5) is a convolution. However, in the most general case 
the POVM does not depend on (0 — 0') only. For ex- 
ample, the conditions (A4) and (A7) are also satisfied by 
all Gaussian POVM's of the type (A8) even when the 
standard deviation a is an arbitrary function of 0'. 



An imperfect detector may be represented by a Positive 
Operator Valued Measure (POVM). To represent imper- 
fect detection of an observable 0, we study the class of 
diagonal POVM's 

n = y #'IW) I 4>'){<t>' I, (Ai) 

where | 0) are eigenstates of 0. These states are assumed 
to constitute a complete set, 

/ #10X01=1. (A2) 

The diagonal form is assumed because interference be- 
tween different detector states should not occur. We as- 
sume that 11^(0') is a nonnegative function. 



APPENDIX B: DERIVATION OF WEAK VALUES 

Prior to the measurement interaction, we assume that 
the object is in the state p s and the pointer is in an ar- 
bitrary mixed state p a . The total density operator prior 
to the interaction has the product form po — p s <S> p a - 

After the interaction, the total density operator has 
evolved to 

p e = UePoUh (Bl) 

where U e is a unitary evolution operator. Since the 
Hamiltonian does not depend explicitly on time, we can 
write 

(j £ = e -ijHdt = e -ieu®P (B2) 
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We may expand this evolution operator to first order in 
e as 



U e w 1 - ie v <g> P. 



(B3) 



To first order in e, the density operator after the interac- 
tion then can be written as 



h = po + ie \po,0®P]- 



(B4) 



After the interaction, we assume that imperfect measure- 
ments are made of the pointer observable Q and the ob- 
ject observable <j>. Each observable is represented by a 
diagonal POVM 



fi = / #'n (0')IWI, 



(B5) 



ftq = j dQ 1 U Q (Q')\Q')(Q'\. (B6) 

Both and IIq are assumed to be classical distribu- 
tions. In this way, the detectors are assumed to be in a 
statistical mixture of pure projector states. We also have 
assumed that the observable <f> has a continuous spec- 
trum. There is no loss of generality in this. 

The joint probability density for <j> and Q after the 
interaction is 



p^Q) =Tr(n n Q/ 3 £ ), 



(B7) 



where the trace is taken over both H s and H a . Using the 
expression (B4) for the first order density operator, we 
find that 



Pe {<t>,Q) = T¥(n p s )TV(n Q/ S a )+ieTr {^U Q \p ,v ® P}) 

= Tr(n^)Tr(n Q/ S a )+2e [Tr(n^ s z>)Tr(n Q p a P) - Tr^vp^Tr^P^) 

I 



(B8) 



We require that the current density of the pointer van- 
ishes, 



(Q I (Ppa + PaP) | Q) = 0. 
This implies that 



Tr 



U Q (Pp a +p a P) 



0. 



(B9) 
(BIO) 



Therefore we can write 

Pe(<P,Q) = Tr(n p s )Tr(nQ/5 a ) 



ie Tr 



%{p s u + up s ) Tr(n Q( 6 a P)(Bll) 



The marginal distribution for the object observable <f> af- 
ter the interaction is 



p e (</>) = / dQp e (^Q). 



(B12) 



Due to the vanishing of the current density of the probe, 
it is found that 



Pe{4>) = Tr(n^/5 S 



(B13) 



Therefore, the probability distribution for the postselec- 
tion observable <fi is unaffected by the measurement in- 
teraction. Note that we have considered an arbitrary 
postselection measurement. This means in fact that the 
probability distribution for every possible object observ- 
able is unaffected by the measurement interaction. 

The conditional probability density for the pointer po- 
sition Q given the postselection outcome <j> is defined as 



p(Q I <P) - 



,Q) 



Pe(0) 



We find that 

p(Q | 4>) = Tr(fi Q/ S a ) + 2^e Re(^)Tr(n Q p a P), (B15) 
where 

Tr(n0j>/3 S ) 



v w {4>) 



Tr(n 0/ 3 S ) 



(B16) 



is the weak value of v for an unsharp postselection rep- 
resented by the POVM fl^. Using Eq. (B9) it can be 
shown that 

(Q | PaP | Q) = \^q{Q I Pa, | Q)- (B17) 

Hence we can write 

Tr(fl Q p a P) = f dQ" [ dQ'n Q (Q')(Q"\Q')(Q'\p a P\Q") 



dQ'n Q {Q')(Q'\p a P\Q') 



= r^j dQ'\\ Q {Q')— <■ n 



"2 I d Q 



d 
dQ 
, dn Q (Q' 



dQ' 



(Q'\p a \Q'), (B18) 



where we have assumed that (Q'\p a \Q') vanishes at the 
integration borders. 

The conditional probability density for the pointer po- 
sition Q given the outcome (p of the postselection is 



(B14) Pe{Q\4>)= j dQ'U Q (Q' + eReu w )(Q'\p a \Q'), (B19) 
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where 
U Q (Q' + eRei> w ) = 



1 + eRe(u w ) 



dQ' 



U Q (Q'). (B20) 



We introduce the expectation value of the pointer posi- 
tion conditioned on the postselection outcome, 

E e (Q | <p) = JdQQ p e (Q | <P). (B21) 

By using Eqs. (B19) and (B20) we find that 

E e (Q | 4>) =E (Q | <f>) + eRe(v w ) 
,dn Q (Q>) 



JdQQ J 



dQ' 



dQ 



^(Qi | p a I Q 1 ). (B22) 



A reorganization of terms gives 
E e (Q | 4>) = E (Q 



eRe(^) 



f dQ 1 (Q'|p a |Q')A J dQQn Q (Q'). (B23) 



We assume that the POVM is unbiased, fulfilling condi- 
tion (A7). This, together with the normalization prop- 
erty of the position distribution, implies that 



E e (Q | <j>) = E (Q | <j>) + eRe{v w ). (B24) 



This shows that the pointer observable Q has been trans- 
lated a distance eRe(is w ), and this setup therefore allows 
for a measurement of Re(v w ). The only restriction on 
the auxiliary pointer system is that the current density 
should vanish. 
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